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Abstract 

The method consists in the monitoring of the temper- 
ature variation of a superconductor tape subjected 
to a trapezoidal current pulse. Simulations of the 
experiment were performed using the lumped capac- 
itance model and also the finite element method. 
The former reproduces the sample's reaction for the 
increasing and decreasing parts of the pulse. The 
same can be said for the finite element method but 
the decreasing temperature as expected is not ad- 
equately reproduced. Treating the problem in two 
dimensions rather than one should correct the situ- 
ation. Nethertheless we obtained a thermal conduc- 
tivity consistent with the literature. 

I Introduction 

Thermal properties are a very good indicator of the 
superconductive transition. However, from a prac- 
tical point of view, the variation of the operat- 
ing temperature is usually very small (say around 



77.4 ± 0.5 K) for power applications. Thus the ther- 
mal properties of high temperature superconductors 
are usually supposed constant when the system is un- 
der normal operating conditions. Actually, the main 
sources of temperature variation are the losses which 
can be separated into a hysteretic contribution (AC 
losses) and other losses that occur because of ther- 
mally activated phenomena, losses in the normal ma- 
trix, ... The true losses are the complete losses (sum 
of those two contributions) and cannot be measured 
by conventional electrical measurements. In order to 
reproduce such a phenomenon without knowing the 
explicit dependence of the true losses with respect 
to external parameters (applied field (DC and AC), 
applied current (DC and AC), temperature, ...), we 
will isolate a part of the sample and then generate 
heat internally. We will then try to modelize the sys- 
tem in order to determine to what degree the thermal 
properties are affected. 
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Fi gure 1; Principle of the experiment 

II Description of the method 

The experiment is very simple and its principle is 
represented in Fig^ 

Part of the BSCCO tape is sandwiched between 
the two halves of a polystyrene block. The joints are 
covered with a silicon adhesive so that no liquid ni- 
trogen will reach the protected part of the tape when 
the block is immersed in the cryogenic bath. 

The data reported here were obtained with a 
BSCCO tape from BICC. A CERNOX sensor from 
Lakeshore Cryotronics was placed at the tape cen- 
ter to measure the temperature inside the block. We 
believe that the higher complexity of the experiment 
using a CERNOX sensor instead of a thermocouple is 
compensated by the absolute temperature measure- 
ment read by the CERNOX sensor (compared to rel- 
ative measurements for thermocouples). 

The temperature of the nitrogen bath is not con- 
stant for two major reasons : the stratification of 
the fluid can change the local temperature and the 
liquefaction of ambient oxygen (with a boiling tem- 
perature of 90K) can change the average temperature 
of the bath. This drift cannot be evaluated with only 
one sensor : ideally, one would include another sensor 
outside the block (soldered to the sample) in order to 
monitor the thermal drift. 

In this article, data from only one sensor are used 
but we are currently using three sensors distributed 
along the protected part of sample. Sensors are 
mounted on small sapphire rectangles placed trans- 
versely across the tape so that we take into account 
all the heat contributions from a given cross-section. 
The thermal drift was substracted from the data us- 
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Figure 2: Power pulses with current ramps of 2 A/s and 
0.2A/s 

ing linear interpolation between the initial and final 
temperatures. 

Our experimental data consist of two measure- 
ments, both obtained from the same sample, in the 
same cryostat and with the same maximum current 
(13 A). The shape of the current pulse applied is 
trapezoidal. In one case the slope is 2 A/s and in the 
other 0.2A/s. The resulting power pulses are repre- 
sented in Fig El 

Below, we shall compare our result for the thermal 
conductivity with those of ref [3] and @] . These de- 
pends strongly on the type of sheath that envelops 
the superconductor. 

The specific heat varies from sample to sample but, 
according to 0, we have, for typical BSCCO tapes, 
0.150 < CpiJ.g^.K- 1 ) < 0.190. A 10% margin is to 
be considered for the same sample because of eventual 
inhomogeneities . 

Ill Numerical Calculations 

In this section, we will address the problem by two 
different methods : firstly the so called "lumped ca- 
pacitance model" and then a more numerical ap- 
proach : the explicit ID finite element method. 

III.l The lumped capacitance model 

For a review of heat transfer, see ref [Q. The main 
hypothesis of this model is that the temperature 
is spatially uniform. This implies that tempera- 
ture gradients within the solid (in this case the 



superconductor) are negligible. Thus we cannot 
use the standard framework of heat conduction 
expressed by the Fourier law (Eq. QJ. Instead, we 
calculate the transient response of the temperature 
T which is determined by formulating the overall 
energy balance of the solid. So if the solid is in 
a fluid (namely nitrogen) with a convection factor 
h {W.mT 2 -K^ 1 ) at a temperature Too, then the 
equation E stora ge = —E ou t becomes : 

jrri 

PVC V — = -hA.iT-T*,) (1) 

A s is the surface of the sample in contact with 
the convective fluid, p (kg.m~ 3 ), V (m 3 ), C p 
(J.kg -1 .K^ 1 ) are the density, the volume and the 
specific heat of the sample respectively and t (s) the 
time. 

In this case, we can define a thermal time constant 
for the system : 

r t =(^)ipV)=R t C t (2) 

where Rt = (o - ) * s ^ ne resistance to convection and 
Ct — pV is the lumped thermal capacitance. 

In order to know if such a regime is reached, 
we must examine the dimensionless Biot number, 
Bi = (hL c ) I (k) where L c is the characteristic length 
of the problem and k is the thermal conductivity 
(W.m^.K -1 ). In order to apply the lumped capac- 
itance model, the Biot number must be very small 
compared to unity. If this is the case, we can say 
that the resistance to conduction within the solid is 
much less than the resistance to convection across the 
fluid boundary layer. 

In our case, the studied system is more complex 
but if the same representation can be used, we can 
modelize our system by the sample, the sample holder 
and the global environment (h, Too). If there are heat 
leaks in the sample holder this will lead to two time 
constants rather than one : one time constant is due 
to the unidimensional conduction (exchange of heat 
at the extremities of the sample) and the other one 
due to the heat exchange between the sample and the 
sample-holder. 

The net outcome of this method is an exponential 
behavior. By fitting the data with exponentials we 
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Figure 3: Fit with two exponentials (lumped capacitance 
model), following the application of the 2A/s trapezoidal cur- 
rent pulse. The upper part represents the error between the 
data and the fit 

implicitly assume that this model is correct and we 
modelize the sample holder action by introducing the 
second time constant. Even if the most interesting 
information is surely the sample's contribution, that 
of the sample holder must be studied if we want to 
ensure that our sample is "well" isolated. Moreover, 
we will be able to determine a characteristic time con- 
stant of our system (sample + sample holder). After 
this time, our system can no longer be considered iso- 
lated since the heat flow through the sample holder 
must be taken into account. 

Another aspect of our problem is that the heat is 
generated internally by an electric current and then 
flows outside the system. If the heat generated is im- 
portant (compared to the leaks) and the heating-time 
is short (which means that the heat exchange between 
sample and sample-holder is small compared to inter- 
nally generated heat) we can neglect the second time 
constant for the heating period only. 

The fit with two exponentials, for the 2A/s case is 
represented in Fig. [3] Similar results were obtained 
with the 0.2A/s pulse. 

III. 2 ID explicit finite element mod- 
elization 

For a review of the finite element method, see ref . 
The master equation is Fourier's law : 

kA s <§ + iV = pVC p <§ (3) 
In order to use the finite element method, the space 



3 



is "discretized" in spatial and time intervals, Ax and 
At. Then we can write : 



, , T p .-7* 

kA ,n - 1 m 



Ax 



rpp _*J>P 

• kA * m+ L - + iA a Ax= (4) 

pA s AxC p Tm A ~ Tm 



Isolating the (T^ +1 ) term, we can write the follow- 
ing equations for all the nodes except the first and 
last ones : 
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We can define the dimensionless Fourier number 
Fo = (^yi and then express the preceding equation 
in term of this number : 



(Ax) 



^+T? n {\-2Fo) = Tl +1 
(6) 



The description for the chain extremities depends 
on the spatial initial and final conditions. In our case, 
we supposed that the sample is in a liquid at a tem- 
perature Too with a convection factor h. After defin- 
ing the dimensionless Biot number (Bi — ^j^-) we 
can write : 



T£ +1 = 2Fo + Bi Too + q^f) + 

(1 - 2Fo-2BiFo)T^ 
and 



(7) 



= 2Fo 
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y 2k 



(1 - 2Fo - 2BiFo)TP 



If we want this equation to have a physical mean- 
ing : (1 — 2Fo — 2BiFo) > 0, we must impose a 
condition on Ax and At. 

Since the discretization of space and time is not 
entirely free, we had to choose our steps carefully. 
The preceding inequality can be expressed in terms 
of Ax and At : 

{Ax) 2 pC p 



At < 



2(k + hAx) 



Of course, this equation can always be fulfilled (by 
diminishing At) but the corresponding increase in 
computing time can become unrealistic. 



Figure 4: Fit with finite element analysis following the ap- 
plication of the 2A/s trapezoidal current pulse 



We implemented this algorithm in C++ (source 
avalaible upon request). 

Fig 0] represents the results of the finite element 
method for the power pulse with a 2A/s current slope 
with the experimental data. Again, similar results 
are obtained with the 0.2A/s pulse. The proposed 
modelization is accurate only for the increasing tem- 
peratures. This is an awaited result because our mod- 
elization is, for the moment, stricly unidimensional. 
As the preceding modelization shows us, there are 
heat leaks in our system and we must take them into 
account. 

Indeed, with our finite element simulation, we take 
into account only one contribution : the thermal con- 
duction of heat along the sample. In order to take 
into account this phenomenon, we have to consider 
the two dimensions of the sample and introduce an- 
other parameter (h!) that will describe the heat con- 
duction along the y axis. 

The simulations give us a thermal conductivity be- 
tween 100 and 125 W.m _1 .K _1 . These results are con- 
sistent with the data cited in references [Hj and 0] 
in which the thermal conductivity is seen to decrease 
with increasing gold concentration. The results for 
our 6% gold sample falls between the 2.9% and 11% 
samples of ref [Hj and between the 4% and 10% sam- 
ples of ref 0]. 
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IV Conclusion 



We have evaluated the thermal conductivity of high 
temperature superconductors with a very simple ap- 
paratus. Whereas the precision is not high (10%), 
it is coherent with other thermal conductivity mea- 
surements made on BSCCO tapes. The use of three 
CERNOX sensors will gives us a better precision and 
will negate the influence of the thermal drift of the 
liquid nitrogen bath. 
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